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13  ABSTRACT 


Thih  report  summarizes  technical  details  of  work  performed  at  the  Ohio  State 
University  ElectroScience  Laboratory  during  the  period  January  8,  1971  to  July  8, 
1971. 

Progress  is  reported  on  the  design,  construction,  and  testing  of  apparatus 
and  computer  software  for  an  experiment  on  atmospheric  imaging  and  image  reconstruc¬ 
tion  at  10.6  microns.  The  initial  experiment  involves  recording,  in  digital  format, 
the  atmospherically  degraded  images  of  a  pair  of  point  sources  and  is  to  be  applied 
later  to  images  received  from  laser  illuminated  real  objects.  Both  experimental 
equipment  data  processing  for  image  reconstruction  are  discussed.  Preliminary  data 
obtained  by  imaging  an  atmospherically  degraded  laser  beam  during  recent  tests  at 
the  RADC  Verona  PATS  test  site  are  presented. 

An  analysis  of  errors  resulting  from  tape  drop-out  during  digital  recording 
using  a  serial  (Miller  Code)  technique  is  presented  for  both  linear  and  logarithmic 
processing  cf  the  analog  image  data  prior  to  recording.  It  is  shovn  that  linear 
processing  is  generally  superior  for  the  present  digital  system  and  hence  is 
recommended. 
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ABSTRACT 


This  report  summarizes  technical  details  of  work  performed  at 
The  Ohio  State  University  ElectroScience  Laboratory  during  the  period 
January  8,  1971  to  July  8,  1971. 

Progress  is  reported  on  the  design,  construction,  and  testing  of 
apparatus  and  computer  software  for  an  experiment  on  atmospheric 
imaging  and  image  reconstruction  at  10.6  microns.  The  initial  experi¬ 
ment  involves  recording,  in  digital  format,  the  atmospherically  de¬ 
graded  images  of  a  pair  of  point  sources  and  is  to  be  applied  later 
to  images  received  from  laser  illuminated  real  objects.  Both  experi¬ 
mental  equipment  and  data  processing  for  image  reconstruction  are 
discussed.  Preliminary  data  obtained  by  imaging  an  atmospherically 
degraded  laser  beam  during  recent  tests  at  the  RADC  Verona  PATS  test 
site  are  presented. 

An  analysis  of  errors  resulting  from  tape  drop-out  during  digital 
recording  using  a  serial  (Miller  Code)  technique  is  presented  for  both 
linear  and  logarithmic  processing  of  the  analog  image  data  prior  to 
recording.  It  is  shown  that  linear  processing  is  generally  superior 
for  the  present  digital  system  and  hence  is  recommended. 
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I.  INTRODUCTION 


This  is  the  fourth  semiannual  technical  report  on  Contract  Number 
F  30602-70-C-0003  entitled  "Investigation  of  10.6  Micron  Propagation 
Phenomena".  This  report  covers  the  period  January  8,  1971  to  July  8, 
1971. 


The  initial  aim  of  this  program  was  to  provide  theoretical  backup 
information  to  the  RADC  Atmospheric  Optical  Propagation  Studies  Program. 
More  recently  this  project  has  centered  on  the  study  of  atmospheric 
imaging  and  restoration  of  atmospherically  degraded  images.  The 
general  aim  is  to  provide  initial  attempts  at  restoration,  to  determine 
potential  problems;  and  to  obtain  initial  data  useful  for  preliminary 
systems  design.  Towards  this  end  a  system  is  near  completion  which 
will  record  10.6  micron  images  and  similarly  the  software  to  restore 
and  examine  these  images. 

Primary  effort  during  the  past  six  months  has  been  devoted  to  the 
design,  construction,  and  testing  of  the  10.6  micron  imaging  system  and 
in  writing  computer  software  for  processing  the  experimental  data. 
Detailed  descriptions  of  the  basic  experiment!  and  of  the  experimental 
apparatus?  have  been  given  in  previous  technical  reports.  Tests  of 
Lhe  analog  portions  of  the  experimental  apparatus  recently  performed 
at  the  RADC  Verona  PATS  test  site  are  described  and  preliminary  results 
of  computer  analysis  of  the  data  thus  obtained  are  presented  in  the 
following  sections.  Changes  to  the  imaging  system  following  the 
Verona  tests  are  also  described  and  an  analysis  of  the  relative  merits 
of  linear  vs  logarithmic  pre-processing  of  the  data  prior  to  digital 
recording  is  included  as  an  Appendix. 


II.  IMAGING  AND  IMAGE  RESTORATION  APPARATUS 

A.  Description 

The  basic  experiment  to  be  performed  following  completion  of  the 
experimental  apparatus  involves  simultaneous  recording  of  two  atmos¬ 
pherically  degraded  laser  beams.!  The  image  received  from  one  beam 
will  be  processed  to  obtain  the  atmospheric  MTF  which  will  then  be 
applied  to  restore  the  degraded  image  received  from  the  second  beam. 
Information  on  isoplanatic  patch  size  will  be  obtained  by  varying 
the  spacing  between  the  two  beams  and  noting  the  effects  on  image 
restoration  and  by  studying  the  correlation  between  the  two  images 
as  a  function  of  beam  separation.  Later  experiments  will  involve 
image  restoration  of  reflected  imaqes  from  real  objects  using  a 
point  source  reference  (laser  beam  or  corner  reflector)  to  obtain 
the  atmospheric  MTF,  and  finally,  the  use  of  edge  effects  in  an  image 
to  eliminate  the  need  for  a  reference  beam  is  to  be  investigated. 
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The  experimental  equipment  described  fully  in  a  previous  report 
consists  of  a  two-channel  scanning-disc  imaging  detector  capable  of 
simultaneous  reception  of  two  separate  images  using  a  40  x  40  element 
array  at  200  frames  per  second  together  with  real  time  digitizing, 
recording,  playback,  and  viewing  equipment.  Single-channel  analog 
portions  of  this  equipment  have  recently  been  completed  and  suc¬ 
cessfully  tested  at  the  RADC  Verona  PATS  test  site  as  described  in 
the  following  section.  The  equipment  has  now  been  returned  to  the 
OSU  ElectroScience  Laboratory  where  the  digitizing  and  recording 
equipment,  the  second  IR  detector,  and  remaining  controls  are  being 
incorporated. 

One  possible  change  iri  equipment  with  respect  to  that  described 
previously?  is  that  the  proposed  logarithmic  video  amplifier  may  be 
omitted  from  the  final  design.  This  device  was  included  in  the  original 
design  when  analog  recording  was  being  considered  in  order  to  improve 
the  signal -to-noise  ratio  of  low-level  signals.  Since  digital  (Miller 
Code)  recording  has  been  implemented,  use  of  a  linear  video  amplifier 
is  found  to  be  generally  superior  (see  Appendix)  as  well  as  simpler  and 
more  reliable  and  hence  is  recommended. 

A  second  change  already  incorporated  into  the  equipment  is  a 
second  set  of  display  electronics,  which,  together  with  a  second  display 
scope  to  be  supplied  by  RADC,  will  permit  simultaneous  viewing  of  both 
video  channels  simultaneously  as  well  as  facilitate  system  alignment 
and  trouble-shooting. 

B.  Preliminary  System  Tests 

Single  channel  analog  video  tests  of  the  scanning  detector  were 
recently  conducted  at  the  RADC  Verona  PATS  test  site  to  verify  major 
system  design  and  performance  criterion  and  to  obtain  preliminary 
data  for  checking  computer  software.  The  equipment  arrangement  used 
for  these  tests  is  shown  in  Fig.  1. 

A  10.6  micron  laser  beam  after  propagating  through  the  turbulent 
atmosphere  is  received  by  the  16  in.  Boller-Chivens  telescope  and 
refocussed  onto  the  scanning  disc  via  the  positioning  and  magnifying 
optics  which  provide  magnifications  of  1,  2,  or  4  as  desired.  Energy 
passing  through  the  scanning  disc  is  converted  to  an  electrical  signal 
by  the  Cd  Hg  Te  detector,  amplified  by  a  preamplifier,  further  amplified 
and  combined  with  blanking  signals  via  a  video  amplifier  and  is  then 
applied  to  intensity  modulate  the  display  and  recording  oscilloscopes. 
Sync  signals  are  also  generated  by  the  scanning  disc  and  sync  generators 
and  are  applied  to  the  raster  generator  which  produces  horizontal  and 
vertical  digital  (point-by-point  staircase)  deflection  voltages  for 
the  video  display  and  recording  scopes  arid  a  blanking  signal  to  turn 
off  the  trace  between  image  elements  and  between  frames.  The  video 
display  scope  provides  a  visual  intensity  modulated  x-y  view  of  the 
infrared  image  received  by  the  scanning  detector,  while  the  vid°o 
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recordinq  scope  provides  a  line-by-line  display  which  is  photographed 
by  a  continuous  motion  strip  camera  for  subsequent  analysis  using  the 
OSU  16  mm.  film  digitizer  as  described  later  in  Section  III -D  of  this 
report. 


Fig.  1.  Experimental  arrangement  for  preliminary  tests 
of  the  scanning  detector. 


III.  DATA  PROCESSING 

The  overall  goal  of  this  program  is  to  study  the  atmospheric 
imaging  and  image  restoration  processes  at  10.6  microns.  The  imaging 
equipment  described  in  Section  II-A  has  therefore  been  developed  to 
record  such  images.  In  Sections  III-C  and  D  the  accomplishments  made 
in  restoring  these  degraded  images  are  described.  First  though,  in 
Sections  I II-A  and  B,  basic  principles  are  reviewed  from  which  the 
equations  used  by  the  computer  software  to  restore  a  degraded  image 
are  derived.  In  Section  III-C  the  computer  software  developed  at 
the  ElectroScience  Laboratory  is  described  and  in  Section  III-D  some 
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preliminary  test  results  obtained  using  this  software  are  presented. 

A.  Image  Restoration 

This  discussion  of  image  restoration  is  limited  to  an  optical 
system  for  which  the  objects  being  imaged  are  incoherently  illuminated 
and  the  point  spread  function  (PSF)  associated  with  each  point  con¬ 
tained  in  a  finite  region  of  the  object  plane  is  assumed  invariant. 
That  is,  the  object  being  imaged  is  assumed  to  lie  entirely  within 
one  isoplanatic  patch.  For  present  purposes  a  single  10.6y  laser 
generated  point  source  is  considered.  Thus,  this  point  source  may 
De  assumed  to  be  either  a  coherent  or  an  incoherent  object  Under 
these  restrictions,  each  point  of  the  image  is  a  weighted  linear 
superposition  of  the  power  distribution  associated  with  that  point 
and  all  neighboring  pcints  in  the  image  plane.  For  the  case  con¬ 
sidered  here,  this  reduces  to  the  power  distribution  resulting  from 
a  single  point  in  the  object  plane.  The  image  irradiance  (watts/m?) 
map  obtained  by  such  a  system  is  described  by  the  convolution  integral 


(1) 


Mx,y)  = 

1  J  J 


I(x'.y')  P(x-x',y-y')  dx'  dy' 


where  I(x',y')  is  the  irradiance  map  of  the  ideal  image,  P(x-x',y-y') 
is  the  PSF  of  the  optical  system  and  Ii(x,y)  is  the  resulting  ir¬ 
radiance  map  in  the  image  plane.  For  an  imaging  system  where  the  area 
between  the  object  and  image  planes  is  void  of  atmospheric  turbulence 
effects,  the  integral  in  Eq.  (1)  represents  the  convolutions  of  the 
telescope  pupil  function  with  the  ideal  image  irradiance  map.  How¬ 
ever,  when  atmospheric  turbulence  is  present  the  system  PSF  becomes  a 
convolution  of  the  pupil  function  with  the  atmospheric  PSF  which  is  the 
case  of  present  interest. 


Let  Ps(x-x' ,y-y' )  be  the  overall  system  PSF  with  atmospheric 
turbulence  effects  present,  then  the  image  irradiance  map  in  terms 
of  the  convolution  integral  is  now 


(2) 


Ij(x,y)  = 


oo 

I (x ' ,y ' )  Ps(x-x',y-y')  dx'  dy'. 

win 


The  function  of  interest  here  Is  the  Ideal  image  irradiance  map, 
I(x',y').  A  solution  for  I(x',y')  is  easier  to  obtain  if  the  Fourier 
Transform  of  Eq.  (2)  is  used.  Such  a  transform  leads  to  the  spatial 
frequency  domain  equation 


(3)  Vw  -  i(Kx.Ky)  Ps(Kx,Ky) 

where  the  indicates  a  Fourier  Transformed  function  and  K  and  K 
are  spatial  frequency  variables.  Here  x  y 
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(4) 


k 


PS(KX,K  )  =  P(VV 

where  A(Kx,Ky)  is  the  frequency  domain  representation  of  the  atmos¬ 
pheric  PSP.  ‘The  spatial  frequency  domain  representation  of  a  PSF  is 
called  the  Optical  Transfer  Function  (OTF).  The  Modulation  Transfer 
Function  (MTF)  is  the  magnitude  of  the  OTF  normalized  by  the  DC 
spatial  frequency  component.  The  MTF  may  be  interpreted  as  giving 
an  indication  of  the  smallest  detail  on  the  object  that  may  be  resolved 
by  the  optical  system. 


Equation  (3)  suggests  that  the  ideal  image  may  be  obtained  by 
dividing  the  transform  of  the  image  irradiance  map  by  the  system 
OTF  and  then  performing  an  inverse  Fourier  Transform 


(5) 


I(x.y)  =  F_1  {I(Kx,Ky)>  =  F-'-* 


Ps(Kx,Ky) 


However,  any  physically  realizable  system  that  might  be  used  to  record 
Il(x,y)  and  Ps(x,y)  will  inherently  introduce  noise  into  these  signals. 
With  the  presence  of  this  noise,  it  is  not  possible  to  obtain  the  exact 
solution  as  indicated  by  Eq.  (4).  There  will  be  components  of  P$(Kx»Ky) 
which  are  zero  or  at  least  very  small  while  the  amplitude  of  these  com¬ 
ponents  contained  in  the  recorded  image  spectrum  will  be  higher  because 
of  noise  components  present  at  these  points.  Calculating  components  of 
the  restored  spectrum  using  Eq.  (4)  at  such  points  will  greatly  amplify 
the  noise  components  resulting  in  an  unintelligible  restored  image. 


Since  an  ideally  restored  image  can  not  be  obtained  in  the  presence 
of  noise,  an  approximately  restored  image  is  obtained  by  replacing  the 
reciprocal  of  the  OTF,  found  from  the  recorded  data,  by  a  processing 
function.  In  the  spatial  frequency  domain,  this  processing  function  is 
equal  to  the  reciprocal  of  the  OTF  spectrum  for  components  which  are 
much  greater  than  the  noise  level,  but  is  such  that  it  greatly  at¬ 
tenuates  the  large  amplitude  noise  components  discussed  above. 

Replacing  l/P$(Kx»Ky)  in  Eq.  (4)  by  the  processing  function  Pp(Kx,Ky), 
and  realizing  £hat  the  resultant  function  will  be  an  approximately  re¬ 
stored  image,  Ia(Kv,K  )  we  have 
m  x  y 


(6) 


‘a<kx>V-  Vw  fyw- 


Further  if  we  assume  the  form  of  P  (K  ,K„)  to  be 

p  x  y 


(7) 


Pp(Kx.Ky)  -  5T-1 


—  P„(K„.KV) 

Ps(Kx.Ky)  y 
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then  by  using  this  processing  function  the  relationship  between  the 
degraded  image,  the  ideal  restored  image,  and  the  approximately 
restored  image  is  given  by 


■v  *  LflC.KJ  p0(kx-k„) 

(8)  IA(Kx,Ky)  =  I(Kx,Ky)  P0(Kx.Ky)  =  1  \  *  0  X  * 


Ps(Kx,Ky) 


-  'i<Kx-V  WV' 


We  can  see  from  Eq.  (8)  that  the  approximate  solution  now  has  resolution 
less  than  that  of  the  ideal  image. 


Taking  the  inverse  Fourier  Transform  of  Eq.  (8)  we  find  the  equation 
to  be  solved  for  the  approximate  restored  image  is 


(9)  IA(x,y)  =  F-1  Uj(Kx,Ky)  Pp(Kx,Ky)}  . 


The  software  developed  at  the  ElectroScience  Laboratory  to  restore 
turbulence  degraded  images  computes  each  component  of  the  approximate 
restored  image  spectrum  using 


(10)  I(KX,K  )  -  Ijd^.K  )/Ps(Kx,K  ) 


(the  inverse  of  Eq.  (5))  and  then  checks  these  values  against  a  test 
value.  This  test  value  is  determined  by  looking  ac  the  entire  spectrum 
found  from  Eq.  (10)  and  choosing  a  level  below  which  all  but  a  few 
obviously  incorrect  components  are  found.  All  components  computed, 
using  Eq.  (10),  which  are  less  than  this  test  value  are  included, 
unaltered,  in  the  spectrum  lA(Kx,Ky);  however,  if  a  component  exceeds 
the  test  value  it  is  considered  to  be  the  result  of  noise  amplification 
as  discussed  previously  and  is  replaced  by  a  more  realistic  value  found 
by  interpolation  of  the  nearby  acceptable  values  of  Ifl(Kx»Ky)*  Thus 
the  P0(Kx,Ky)  factor  of  the  processing  function  (Eq.  (7))  is  equal  to 
1.0  when  the  spectrum  components  are  less  than  the  preselected  test 
value  and  equal  to  some  complex  constant  a,  where  |a|  <<  1.0,  for 
those  components  that. exceed  the  test  value.  That  is 


(ID 


P0(Kx>V  ■ 


1 .0  for 


lMy> 
P'sOvy 


<  Test  value 


ct(  | a  |  <<  1 )  for 


Vvv 


>  Test  value. 
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Substituting  Eq.  (11)  into  Eq.  (7)  we  find  the  processing  function 
to  be 

\ -  for  spectrum  components  <  test  value 

PS<Kx’Ky) 


-  for  spectrum  components^  test  value. 

Ps(Kx,Ky) 

This  function  does  indeed  behave  according  to  the  format  previously 
stated  as  desirable  for  the  processing  function. 

We  are  thus  interested  in  solving  Eq.  (9)  where  Pp(Kx,Ky)  is 
given  by  Eq.  (12) 


B.  The  Fast  Fourier  Transform 


The  method  of  solution  proposed  in  Eq.  (9)  of  the  preceding 
section  indicates  that  a  capability  for  performing  a  bidirectional 
Fourier  Transform  is  required.  Since  the  ElectroScience  Laboratory 
imaging  hardware  generates  1600  discrete  samples  per  image  (40  lines/ 
image  and  40  samples/line)  the  transform  that  must  be  applied  is  the 
(finite)  Discrete  Fourier  Transform  (DFT).  The  defining  equation 


for  a  two  dimensional  DFT  is 

N  -1  N  -1 

(13)  >u  =  fnr  Z  ? 

Vy  R=0  S=0 


For  the  ElectroScience  Laboratory  system  N  =  N  =  N  =  40  samples. 

X  y 

The  validity  of  this  transform  pair  can  be  proven  by  substituting 
1 1 j  from  Eq.  (13)  into  Eq.  (14)  and  using  the  orthogonality  relation¬ 
ship  for  exponentials 
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V1  V1 

(15)  l  l  e 

1=0  J=0 


-«■(?  *f)  »■(?•  f)_ 


N  N  for  R=M 
x  y  S=N 


0  otherwise. 

Such  a  DFT  pair  will  be  denoted  as 
If 


(16) 


RS 


IJ  * 


Two  other  properties  of  the  DFT  pair  worth  mentioning  here  are  first 
that  Irs  and  Iy  are  periodic  with  periods  of  ±kNx  and  ±s%  (k, £=0,1 ,...), 
and  second,  that  the  transforms  of  two  real  valued  sequences  are  obtained 
simultaneously  by  forming  a  sample  by  sample  complex  sum  of  the  two  real 
sequences  and  then  applying  the  DFT  to  this  complex  sequence. 

Although  the  DFT  is  a  very  powerful  bidirectional  mapping  operator, 
it  is  very  time  consuming  and  hence  it  is  costly  to  use  these  transforms 
directly  when  a  large  number  of  samples  are  involved  as  in  the  case  here. 
To  directly  apply  the  DFT  to  the  1600  data  samples  from  a  single  two 
dimensional  image  would  require  =  2.56  x  106  complex  multi  pi  ication 
and  addition  operations.  Thus,  performing  the  TFT  transform  necessary 
to  find  a  restored  image  becomes  a  natural  application  for  the  Fast 
Fourier  Transform  ( FFT ) .  The  FFT  is  a  highly  efficient  method  of 
iteratively  calculating  the  coefficients  of  the  DR  by  combining 
progressively  larger  weighted  sums  of  the  sampled  data.  The  approxi¬ 
mate  number  of  complex  multiplications  and  additions  necessary  for 
calculating  a  two  dimensional  FFT  is  4N2  log2  N(3).  However,  for  the 
FFT,  the  number  of  samples,  N0,  must  be°an  integral  power  of  two. 

Thus  to  include  all  40  x  40  samples  generated  by  the  ElectroScience 
Laboratory  equipment,  these  40  x  40  samples  must  be  placed  in  the 
center  of  a  64  x  64  sample  input  to  the  FFT  with  samples  not  in  the 
40  x  40  array  being  set  equal  to  zero.  The  number  of  operations  . 

necessary  to  calculate  the  FFT  for  N0  =  64  is  approximately  9.8  x  10, 
or  less  than  4%  of  the  operations  necessary  to  calculate  the  DFT  directly. 
Cochran3  et  al  also  point  out  that  not  only  are  the  number  of  calculations 
and  hence  the  computational  time  reduced  by  this  factor,  but  also  round 
off  errors  are  reduced  by  the  same  factor. 

The  particular  FFT  algorithm  used  here  is  one  supplied  in  IBM's 
F0RTLIB  library;  the  subroutine  is  called  F0UR2  and  is  basea  on  the 
Cooley-Tukey  FFT.  F0UR2  computes  the  DFT  of  a  single  precision  input 
data  array  where,  for  the  present  case  of  interest,  this  array  repre¬ 
sents  samples  taken  from  one  cycle  of  the  periodically  extended  two- 


8 


dimensional  image.  By  merely  changing  one  parameter  in  the  calling 
statement  for  F0UR2  either  the  DFT  or  inverse  DFT  is  computed. 

C.  Computer  Software 

Computer  software  ior  restoring  turbulence  degraded  images  has 
been  developed  at  the  ElectroScience  Laboratory.  This  software  is 
written  in  Fortran  IV  G  level  for  the  Ohio  State  University  IBM 
360/75  computer.  The  program  is  not  yet  finalized  since  there  will 
undoubtedly  be  some  changes  required  once  actual  digital  data  from  the 
ESL  Miller  Code  system  is  available  for  analysis.  However,  the  general 
format  described  here  for  the  computer  processing  steps  should  not 
require  changes. 

The  basic  program  format  is  given  in  Fig.  2.  A  pair  of  sampled 
images  are  first  read  in  where  the  PSF  is  denoted  as  Image  1  and  the 
degraded  image  as  Image  2.  These  40  x  40  discrete  sample  arrays  are 
each  centered  in  a  64  x  64  array  (the  FFT  subroutine  requires  that 
each  array  dimension  equal  an  integral  power  of  two)  and  the  remaining 
elements  are  set  equal  to  zero.  The  center  of  gravity  of  each  image 
is  then  calculated  from  which  two  parameters  representing  the  dif¬ 
ference  in  the  x  and  v  coordinates  of  the  two  centers  of  gravity  are 
obtained 

(17)  TN2  =  shift  in  x-coordi nates 

(18)  TNI  =  shift  in  y-coordi nates. 

These  parameters  are  later  used  to  alter  the  phase  of  the  Image  2 
spectrum  so  that  the  transformed  Image  2  will  have  its  center  of 
gravity  at  the  same  location  as  the  PSF  (Image  1). 

The  complex  sum  of  the  two  real  images  is  formed  and  the  FFT  is 
taken  of  this  rum.  The  resulting  spectrum  is  then  unraveled  to  obtain 
the  spectrum  of  Image  1  (i.e.,  the  OTF)  and  that  of  Image  2.  The 
spectrum  of  the  restored  image  is  next  calculated  by  using  Eq.  (10) 
and  testing  to  see  that  no  component  exceeds  the  test  value  (see 
Section  I I I-C. ) .  When  a  component  does  exceed  this  value,  it  is 
replaced  by  a  more  realistic  value  derived  from  the  valid  neighboring 
lower-frequency  spectrum  components.  The  result  is  an  approximate 
spectrum  for  the  restored  image  as  described  previously. 

The  resulting  spectrum  is  transformed  back  to  the  x-y  spatial 
coordinate  system  and  an  IBM  subroutine  PL0T3D  then  outputs  this 
data  in  the  proper  form  to  be  used  as  input  data  to  the  SDS  910 
plotter.  The  result  is  a  three-dimensional  plot  of  the  restored 
image.  For  visual  comparison  purposes,  the  degraded  image  and  the 
system  PSF  are  also  displayed  similarly  in  three-dimensional  plots. 
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Some  of  the  particular  areas  of  the  program  that  may  be  altered 
at  a  later  time  are  the  read-in  and  display  methods  and  possibly 
the  form  of  the  processing  function.  It  is  also  possible  that  in¬ 
formation  other  than  the  restored  image  and  its  spectrum  may  be 
desired  later.  Quantities  such  as  the  system  MTF  and  the  correlation 
between  restored  images  as  a  function  of  time  or  horizontal  dis¬ 
placement  (for  the  twin  beam  experiment)  are  now  available  from  the 
general  program  when  desired. 

In  the  following  section  of  this  report  the  general  program 
format  is  adapted  to  the  specific  task  of  processing  data  scanned 
and  digitized  from  degraded  images  recorded  on  photographic  film. 

D.  Preliminary  Test  Results 

Preliminary  tests  of  the  ElectroScience  Laboratory  10. 6u  scanning 
detector  were  performed  at  the  RADC  Verona  PATS  test  site  during  April 
and  May  1971  as  described  in  Section  II-B.  Using  the  RADC  1000  foot 
outdoor  propagation  range  and  a  10. 6y  laser  source,  the  equipment  was 
set  up  to  display  and  photograph  the  turbulence  degraded  laser  spot 
(see  Fig.  1).  Only  the  video  intensity  and  horizontal  sweep  voltages 
were  applied  to  the  oscilloscope  for  photographic  recording.  The 
vertical  sweep  was  supplied  by  the  constant  film  motion  traveling 
from  top  to  bottom  with  respect  to  the  display  scope.  Although  the 
human  eye  sees  the  crt  display,  with  the  vertical  sweep  applied, 

(Video  display  scope  of  Fig.  1 )  as  a  40  x  40  element  picture  of  a 
turbulence  degraded  spot,  because  of  the  vertical  interlace  used 
by  the  scanning  disc  and  raster  generator  this  image  is  actually 
generated  by  displaying  the  odd  scan  lines  of  the  image  on  a  vertical 
down-scan  ( i . e. ,  1,  3,  ...»  39)  followed  by  the  even  scan  lines  on  a 
vertical  up-scan  (i.e.,  40,  38,  ...»  2).  Thus  when  viewing  the  data 
recorded  by  the  continuous  motion  camera  we  see  two  consecutive  half¬ 
images,  one  essentially  the  mirror  image  of  the  other,  which  compose 
a  single  full  image. 

Upon  receiving  the  filmed  data  at  the  ElectroScience  Laboratory, 
it  was  reviewed  visually  to  select  several  sections  of  interest  to  be 
used  as  test  data  for  the  image  processing  computer  software.  The 
selected  sections  of  film  were  then  digitized  using  the  ESL's  digital 
film  analyzer  which  scanned  each  half-image  and  recorded  the  resulting 
array  of  6-bit  digital  samples  on  IBM  7-track  magnetic  tape.  This 
tape  then  served  as  the  input  to  a  modified  version  of  the  image 
restoration  computer  software. 

Briefly,  these  software  modifications  consist  of  a  series  of 
subroutines  which  read  the  half-frame  data  off  the  magnetic  tape. 
Interlace  the  two  half-images  to  form  a  single  40  x  40  sample  image, 
and  then  convert  this  photographically  negative  image  to  a  photo¬ 
graphically  positive  image. 
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The  first  test  of  the  computer  software  was  performed  by  using 
the  same  image  as  both  the  degraded  image  and  the  system  PSF.  Such 
a  restoration  physically  means  that  the  point  source  of  unity  mag¬ 
nitude  from  which  the  PSF  was  derived  is  being  restored.  We  would 
expect  this  restored  image  to  be  a  delta  function  of  unity  magni¬ 
tude  and  indeed  this  is  the  result  obtained  using  the  image  processing 
software.  Figure  3  is  a  three-dimensional  intensity  plot  of  this 
restored  point  source.  Also  included  are  the  PSF  and  the  degraded 
image  respectively,  in  Fig.  4  and  Fig.  5.  Note  that  these  are 
actually  the  same  function  in  this  instance. 

The  next  test  was  conducted  by  choosing  a  particular  image  as 
the  PSF  and  then  restoring  a  degraded  image  recorded  1/200  sec  later 
in  time.  An  example  case  is  shown  below.  Figure  6  is  the  system 
PSF,  Fig.  7  is  the  degraded  image,  and  Fig.  8  is  the  resulting 
restored  image  obtained  using  the  ESL  image  processing  software. 

Figures  9a  and  9b  show  the  principle  cuts  through  the  restored 
image  in  Fig.  8.  These  principle  cuts  clearly  show  the  restored  image 
in  the  presence  of  noise. 

Obviously  this  test  case  does  not  represent  the  true  method  of 
image  restoration  since  in  the  actual  situation  the  degraded  image  to 
be  restored  and  the  system  PSF  would  be  recorded  simultaneously. 
However,  the  results  above  do  indeed  present  verification  that  the 
computer  software  is  functioning  properly. 


IV.  SUMMARY 

This  report  has  described  work  performed  during  the  last  six 
months  in  construction  and  testing  of  equipment  and  in  writing  computer 
software  for  experiments  in  the  10.6  micron  imaging  and  image  restoration. 
Tests  recently  completed  at  RADC  of  the  scanning  detector  and  analog 
video  display  are  described,  and  samples  of  the  resultant  data  for  an 
atmospherically  degraded  laser  beam  both  before  and  after  restoration 
using  the  ElectroScience  Laboratory  computer  software  are  presented. 

Work  now  in  progress  for  completing  the  experimental  system  and 
experiments  then  to  be  performed  are  summarized. 


12 


-  -  ; 


13 


Fig.  3.  Restored  unit  intensity  point  source. 
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Principle  x  cut  thru  restored  image  in  Fig. 


Fig.  9b.  Principle  y  cut  thru  restored  image  in  Fig. 


APPENDIX 

ERROR  ANALYSIS  FOR  DIGITAL  RECORDING  USING 
LINEAR  OR  LOGARITHMIC  PREPROCESSING  OF  THE  IMAGE  DATA 


The  relative  merits  of  linear  or  logarithmic  amplfication  of  the 
video  signals  from  the  preamplifier  (see  Fig.  1)  prior  to  recording 
are  considered  in  this  appendix.  Two  types  of  system  errors  are 
derived  and  compared  for  a  digital  (Miller  Code)  recording  system 
using  linear  amplification  and  an  equivalent  system  using  logarithmic 
amplification.  First  a  system  is  considered  where  the  only  error 
results  from  a  bit  being  dropped  during  recording  or  playback.  Next, 
this  same  system  is  considered  with  respect  to  the  aperture  error 
associated  with  analog  to  digital  conversion. 

I.  DIGITAL  SYSTEM  DROP-OUT  ERRORS 

The  Computer  Laboratory  Model  HS  802  8  bit  Analog  to  Digital 
(A/D)  converter  used  in  the  ElectroScience  Laboratory  system  has  a 
maximum  unipolar  input  voltage  range  of  2048  mv.  Thus,  the  weight 
of  a  least  significant  bit  (LSB)  is 

(19)  1  LSB  =  M8™.  =  g  mv. 

28 

Let  W  denote  this  weighting  factor  and  Vfl/n  denote  the  maximum  input 
voltage  to  the  A/D  converter.  A/u 

It  has  been  found  from  previous  tests  on  the  ESL  system  that 
the  minimum  signal,  due  to  black  body  radiation  from  the  scanning 
disc,  at  the  output  of  the  Perry  preamp  (see  Fig.  1)  is  approximately 
0.35  mv  peak  and  the  maximum  signal  output  of  this  preamp  is  1.4  volts 
peak.  Thus  let 

(2°)  vmax  =  1.4, 

and 

(21 )  vm.  =  0.35  mv. 

min 

A.  Digital  System  with  Linear  Amplification  (Fig.  10a) 

Let 


(22)  Vc  -  k  vc 

be  the  correct  output  from  the  Miller  Code  decoder  and  let 
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LINEAR  SCHEME 


v  aLOGv+£  S  =  k0  (a  LOG  v  +/3 ) 


Fig.  10.  Block  diagram  of  linear  and  logarithmic 
amplification  schemes. 
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(23)  Ve  =  Rve 

be  an  incorrect  output  from  the  same  decoder  when  a  bit  has  been 
dropped  from  Vc.  Amplifier  gain  k  is  adjusted  so  that  for  an 
input  of  v  =  Vmax  the  amplifier  output  is  V^p.  Thus 

(24)  k  =  I^/D  . 

max 

Now  consider  an  8-bit  digital  word 


I  MSB  I  I  I  II  I  ILSBI 
Bit  Number  8  7  6  5  4  3  2  1 

Weight  128  64  32  16  8  4  2  1 

Voltage  Value  W27  W26  W25  W24  W23  W22  W21  W2°  . 

If  the  signal  Ve  is  a  digital  output  resulting  from  the  loss  of 
the  bit  with  a  voltage  value  of  W2n  then  the  error  is 

(25)  E]  =  Vc  -  Ve  =  Vc  -  (Vc  -  W2n)  =  W2n. 

In  terms  of  vc  and  ve  this  is 

(26)  ,  V(;  .  Ve  =  ,2  =  W|_ 

where  n  =  0,  1,  2,  ...»  7. 

B.  Digital  System  with  Logarithmic  Amplification  (Fig.  10b) 
Let 

(27)  Sc  =  kQ  (a  log  vc  +  e) 

be  the  correct  output  from  the  Miller  Code  decoder  and  let 

(28)  se  =  k0  (“  lo9  ve  + 


23 


be  an  incorrect  output  from  the  same  Miller  Code  decoder  when  a  bit 
has  been  dropped  from  Sc.  Amplifier  gain  k0  is  adjusted  so  that  for 
an  input  of  v  =  vmax  the  log  amplifier  output  is  V^D.  Thus 


(29) 


VA/D 

d  '°9  vmx  +  8  ‘ 


If  this  signal  Se  is  a  digital  output  resulting  from  the  loss  of  the  bit 
with  a  voltage  value  of  W2n  then  the  error  is 


(30)  E0  -  Sc  -  Se  =  Sc  -  (Sc  -  W2n)  = 

In  terms  of  vc  and  ve  this  is 

(31)  EQ  =  ko(a  log  v£  +  e)  -  kQ(a  log  vg 

(32)  Eq  =  kQa  log  vc/ve  =  W2n 

°r  W2n 


(33) 


v  =  v  10 
e  c 


koa 


and  hence  the  error  in  terms  of  v„  and  v„  is 

c  e 


(34) 


A,  =  -  VQ  =  V„ 

o  c  e  c 


/  M2n\ 

(l  -  10 


+  6)  =  W2n 


where  n  =  0,  1 ,  •••  7. 

C.  Comparison  of  Drop-out  Errors  for 
Linear  and  Logarithmic  Systems 

For  the  digital  system  with  error  resulting  from  the  loss  of  a 
bit  during  decoding  wo  have  for  linear  amplification 


(35) 


A  -  W2 

A!  ‘  TT~ 


and  for  logarithmic  amplification 

-  W2n  \ 


(36) 


■’lO- 


1 0  kna 


Oa  J 
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A  comparison  between  ai  and  Ao  will  now  be  made  where  the  log 
amplifier  chosen  is  described  in  Texas  Instrument’s  Bulletin  No 
DL-S  7111427  (page  6).  This  amplifier  has  an  69  dB  range  and  the 
transfer  characteristic  is  plotted  here  in  Fig.  11.  From  this  curve 


1 

8 


1 

1 


Next  the  gains  k  and  k0  are  calculated  to  be 


we  find 

(37) 

a  = 

(38) 

B  = 

Next  the 

(39) 

k  = 

and 

(40) 

V 

Hence  we 

have 

(41) 

A1  : 

and  for 

the  li 

(42) 

Ao  : 

k- W-l.  46 


max 


■A/D  _  2.048  _  - 

■ww?  -  -TO  - 3-34- 


=  ^  -  -*■  *  6.47  x  10"3  2n 


W2n  \ 

=  -  10  ’  -  vin  (l  -  10'16’2  X  10’3  2") 


These  curves  are  plotted  in  Fig.  12  for  Vin  =  vmax  and  for 
vin  =  0.1  vmax.  Prom  Fiq.  12  it  is  seen  that  for  the  lower  bits 
(i.e.,  n  =  0,  1,  2,  3,  4)  the  error  associated  with  the  logarithmic 
system  is  greater  than  that  for  the  linear  system  until  the  input 
signal  drops  to  approximately  0.1  vmax  or  lower.  Note  that  the 
apparent  improvement  for  the  higher  order  bits  when  using  the  log 
amplifier  for  the  low  input  signal  levels  as  indicated  by  the 
curves  are  not  applicable  when  considering  bit  dropout  since, 
for  low  signal  levels,  the  high  order  bits  are  not  used.  Hence 
the  improvement  obtainable  by  using  the  logarithmic  amplifier,  even 
at  low  signal  levels,  is  slight.  The  linear  amplifier  is  therefore 
recommended  for  the  present  system. 
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12  3  4  5 

n*  BIT  DROPPED 


Fig.  12.  Error  voltage  vs  bit  dropp'd  for 
linear  and  log  systems. 


II. 


APERTURE  ERROR 


Aperture  error  for  an  A/D  converter  is  defined  as 
(43)  E  =  ta  ^ 

where  ta  is  the  maximum  aperture  time  and  v  is  the  input  to  the  A/D 
converter. 

For  the  linear  amplification  scheme  (Tig.  10a),  the  input  to  the 
A/D  converter  is 


(44)  V  =  kv 

and  the  aperture  error  is  therefore 

(45)  41  =  ta  ar  =  ta  k  HT 
assuming  constant  gain. 

Next  consider  the  aperture  error  associated  with  the  logarithmic 
amplification  scheme  (Fig.  10b)  where  the  input  to  the  A/D  converter 
is  given  by 

(46)  S  =  kQ  (a  log  v  +  e). 

Aperture  error  for  this  case  is 

(47)  4  ■  ta  £  ■  k0  t a°lSS-e  dv 

again  assuming  constant  gain.  Now  consider  the  case  for  which  the 
amplifier  input  is 

(48)  v  a  1.4  sin  u)Qt 
where 

(49)  w  -  2ir  x  1.6  x  105. 

This  represents  the  highest  amplitude  and  highest  frequency  component 
expected  as  an  input  in  the  ElectroScience  Laboratory  system. 

The  Computer  Laboratory  Model  HS  802  8-bit  A/D  converter  has  a 
maximum  aperture  time 
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(50)  ta  =  10~9  sec 

thus  for  the  linear  scheme  the  aperture  error  is 

(51)  T  =  1.4  ta  k  cos  wQt 
from  which 

(52)  Hlmax  =  1.4  ta  k  (jjq  =  2.04  mv 

while  for  the  logarithmic  scheme  the  aperture  error  is  given  by 

cos  t*  t 

(53)  $  =  k  to  a  oj  log  e  - r- 

o  o  3  sin  4jQt 

Equation  (53)  indicates  that  the  maximum  value  of  $  can  be  infinite 
if  a  sample  is  taken  at  the  wrong  time,  i.e.,  when  the  input  voltage 
is  passing  through  zero.  This,  of  course  does  not  happen  in  practice 
since  no  realizable  amplifier  remains  logarithmic  for  zero  input 
voltage.  For  a  practical  logarithmic  amplifier,  aperture  error  would 
in  general  be  greater  for  low  level  input  signals  and  somewhat 
smaller  for  high  level  signals  than  when  a  linear  amplifier  is  used. 
Since  aperture  error  for  the  linear  system  is  approximately  1/4  the 
value  of  the  least  significant  bit,  this  error  should  be  negligible 
for  either  the  linear  or  logarithmic  amplifier  in  the  present  system. 
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